Tracking quintessence and /c-essence in a general cosmological background 



Rupam Das, Thomas W. Kephart, and Robert J. Scherrer 
Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235 

We derive conditions for stable tracker solutions for both quintessence and fc-essence in a general 
cosmological background, H 2 oc /(/£>). We find that tracker solutions are possible only when r\ = 
din //din p ~ constant, aside from a few special cases, which are enumerated. Expressions for the 
quintessence or fc-essence equation of state are derived as a function of n and the equation of state 
of the dominant background component. 
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I. INTRODUCTION 

The universe appears to consist of approximately 30% 
nonrelativistic matter, including both baryons and dark 
matter, and 70% dark energy (see Ref. [ij for a recent 
review, and references therein) . The evolution of the dark 
energy density depends on its equation of state, which is 
usually parametrized in the form 



PDE = wpoE, 



(1) 



where pbe and pbe are the pressure and density of the 
dark energy. Then the density of the dark energy scales 
as 



Pde oc 



(2) 



The simplest model for the dark energy is a cosmolog- 
ical constant, for which w — — 1 and pde = constant. 
More complex models have been proposed, in which the 
dark energy arises from a scalar field c/>; these are called 
quintessence models 0, IE 0, IE E3 These models gen- 
erally give rise to a time-varying w^, and more complex 
behavior for pde- One advantage of such models is that 
certain classes of quintessence potentials lead to tracker 
behavior, in which the evolution of the scalar field is in- 
dependent of the initial conditions. The conditions for 
such tracking behavior have been worked out in detail by 
Steinhardt, et al. @|. 

A second class of models generalizes quintessence to 
allow for a non-standard kinetic term. These models, 
dubbed fc-essence, have also been explored in great detail 

BIlUHmillHEim- These models can also 
lead to tracking behavior, and the conditions necessary 
for such behavior have been discussed by Chiba [T3 |. 

Both quintessence and fc-essence can be generalized to 
modified versions of the Friedmann equation. In the stan- 
dard Friedmann equation, the relation between the scale 
factor a (or, alternatively, the Hubble parameter H ) and 
the density is 



H 2 = 



P 
3' 



(3) 



where we set 8irG = 1 throughout. However, various 
proposals have been put forward to modify this equation 
at high energy. In type II Randall-Sundrum models, for 
example, one has [lfiL fl7l| 

H 2 cx p\ (4) 



in the limit of large p, while Gauss-Bonnet models can 
give 



H 2 cx p 2 l\ 



(5) 



The Cardassian model 
the form 



with n < 2/3. 



H9| assumes an expansion law of 



Bp" 



(6) 



Motivated by these examples, numerous authors have 
examined the evolution of various dark energ y m odels in 

the context of non-standard expansion laws he m m 

HE l24l l2p . The most general treatments are given in 
Refs. |23| -|25j- Sami et al. examine quintessence 
with constant for a power-law modification to the 
Friedmann equation (H 2 cx p q ). Copeland et al. [24[ dis- 
cuss "scaling" quintessence models, i.e., models for which 
Wcl> = Wbi with an arbitrary expansion law, H 2 oc /(p). 
Here wb is the ratio of pressure to density for the dom- 
inant, "background" fluid, e.g., wb = for a matter- 
dominated universe, and wb = 1/3 for a radiation- 
dominated universe. Tsujikawa and Sami |25| examine 
arbitrary scalar field models (including both quintessence 
and fc-essence) with scaling behavior [w^ = wb) in mod- 
els with a power-law modification to the Friedmann equa- 
tion, H 2 oc p q . 

Here we generalize this earlier work by examining 
tracking solutions for both quintessence and fc-essence 
in a general cosmological background characterized by 
H 2 oc f(p). Although we adopt the approach of Stein- 
hardt et al. [|| for quintessence and Chiba 12] for fc- 
essence, our formalism encompasses tracking solutions 
not only for a wide range of potentials but also for a 
wide range of f{p). We derive sufficient conditions for 
both V(c/>) and /(/?) to obtain tracking solutions with a 
constant w^. This formalism provides us with a generic 
method to study these solutions for a wide variety of 
scalar field models such as quintessence, tachyon, fc- 
essence, and phantom models. 
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II. QUINTESSENCE 
A. Tracking solutions 

The equation of motion for the ^-field is 

4> + ZH4> + v = o, 

where 



V = dV/d<p, 



and 



H 2 



/(/>)• 



(7) 



(8) 



(9) 



Here a is the Robertson- Walker scale factor, and p is the 
total density, given by 



P = Pb 



(10) 



where ps is the background (radiation + matter) density, 
and prf, is the scalar field energy density. The standard 
Hubble expansion law corresponds to equation with 
f(p) = P'i m this paper we allow f(p) to have an arbitrary 
functional form. 

By definition, the tracking solutions are the solutions 
to which the evolution of the scalar field 4> converges for 
a wide range of initial conditions for <j> and <j>. We fol- 
low the approach prescribed by Steinhardt et al. 6] for 
quintessence, but now generalize it to the arbitrary ex- 
pansion law given by equation @ ■ For tracking solutions, 
w<p is nearly constant @|, where w^, is given by 



w$ = — 



I 

2 



<$> 2 -V 



V 



It follows from equation @ that 

H = ^H 2 V [(w^ - w B )[l - fy) - (1 + W4,)], 



(11) 



(12) 



where r\ encodes the information on the generalized ex- 
pansion law in equation ©: 



V = 



dlnf(p) 
ffln(p) ' 



(13) 



For the standard Hubble expansion, 77 — 1. In this pa- 
per, we will confine our attention to the case 77 > 0, and 
our conclusions will be valid only for this case. How- 
ever, we note that 77 < can lead to interesting types 
of behavior (e.g., a phantom- like future singularity in a 
matter-dominated universe). Several specific models of 
this type are mentioned in Ref. po| . 

By combining these relations, it is useful to cast the 
equation of motion into the following form: 



V 



±3H\ 



1 



^(1 + -) 
2 1 6 ' 



(14) 



where x = (1 + w^) / (1 — w^) = \§ 2 jV is the ratio of the 
kinetic to potential energy for </>, and x' = d hix/d In a. 
The ± sign depends on whether > or V^, < 0, re- 
spectively. It follows from equation (|14fl that the tracker 
condition (itJ^ w 0) becomes 



V 



(15) 



This is the generalization of the Steinhardt et al. [6( 
tracking condition to an arbitrary expansion law. 
As in Ref. we define the function 

T v = V^V/iV^f, (16) 

whose properties determine whether tracking solutions 
exist. By taking the time derivative of equation l|14l) 
and combining with the equation (|12fl and (|14|) itself, we 
obtain the following equation: 

-p 1 + V _ r)(w B - w^VLb 



2 2(1 + ^) 

t](wb - w^Qb + T) + (tj - 2)w(jj x' 



2(1 + w ) 



Q + x' 



(I + W4) (6 + x') 



!\2 ■ 



(17) 



where x" = d 2 Inx/d In a 2 . As expected, equation 117() 
reduces to the corresponding equation in Ref. for 
77 = 1 . In a universe dominated by a background fluid 
(Qb ~ 1) with tx>0 w constant and nearly constant IV, 
the above equation becomes 

_ V + 1 v(w B - w^) 
V ~ 2 + 2(1 + ^) ' 

1 77 ( 1 + Wb 

2 + 2 



1 



Wcj, 



(18) 



In deriving the above equation, the plausibility of the 
condition that Ty ~ constant has been discussed in de- 
tail in Ref. |6j. The crucial point is that this condition 
encompasses a wide range of potentials including inverse 
power law potentials and combinations of inverse power 
law terms to give rise to tracking solutions. 

We must know the appropriate restrictions on 77, i.e., 
on /(/?) to extract the tracking solutions from equation 
(|T7jl . Since the left-hand side of equation ((TH1) is nearly 
constant, it follows that 77 must be nearly constant during 
background domination, i.e., the function /(/?) must sat- 
isfy (|13fl for a nearly constant 77. Thus we require an extra 
condition, in addition to the conditions on Ty, to derive 
tracking solutions for both quintessence and fc-essence. It 
is obvious that this extra condition arises from the extra 
"degree of freedom" in choosing a different cosmological 
background. The only case for which 77 is exactly constant 
is f(p) oc p n for a constant n. This power-law behavior 
includes both the Randall- Sundrum and Gauss-Bonnet 
models as special cases, and it was studied in detail in 
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Ref. [23j. Of course, more general conditions can pro- 
duce an expression for f(p) that is roughly constant over 
a wide range in the scale factor. For instance, a sum of 
power laws, e.g., as in equation ©, gives a value for 77 
that is nearly constant over most of the evolution of the 
universe, i.e., at all times except for the epoch when the 
two contributions to f(p) are roughly equal. 

Note that there are a few trivial special cases for which 
this argument breaks down. In particular, if V is a 
constant, the right hand side of equation (|Tlf> must be 
zero; this can be achieved by taking w<p = ±1. The case 
w<f, = — 1 corresponds to a non-zero constant potential, 
while Wtf, — 1 is the solution for V — 0. Both of these re- 
sults are independent of the value of H on the right-hand 
side of equation Ijl4|l and are therefore independent of 77. 

The validity of equation (|18|) may be checked by com- 
paring with the results obtained by Sami et al. [23| . 
For scaling solutions with a constant w$ in a background 
dominated universe, the potential function takes the fol- 
lowing form [2^ 

V{(j)) cx (T Q , (19) 
where a is constant. Then we obtain from equation 118|) 



1 + W0 « 77(1 + W B ) 



a + 2 



(20) 



This solution agrees with the result obtained in Ref. [2c 



B. Stability of the tracking solutions 

So far, we have derived solutions with constant in a 
general cosmological background; now we want to check 
the stability of these solutions with constant uu. In order 
to check the stability, we perturb the tracker value of w<f, , 
which we will call wq, by an amount 6. Then we expand 
equation (|17fl to lowest order in S and its derivatives to 
obtain 



25 +3[j]{l + w B )-2w ]S 
+9t 1 (1 + w b )(1-w )S = Q, 



(21) 



where the prime means d/dhia and wq is the value of 
derived from equation (jT%|l . The solution of this equation 



(22) 



where 



7 = -- [77(1 + w B ) - 2w ] 



±^877(1 + w B )(l - w ) - [r?(l + w B ) - 2w } 2 . 

(23) 

In the derivation of this equation, Ty and 77 are assumed 
to be constant. 



In order to have 5 decay, the real part of 7 has to be 
negative. Hence, it follows that 



w < 



7/(1 + W B ) 



(24) 



provided the quantity under the square root is positive. 
If the quantity under the square root is negative (so that 
both values are real), then the above equation is also a 
necessary condition since the first term under the square 
root is always positive, provided 77 > and wq < 1. 
Using equation l(TTj|l , the above inequality can be written 
in terms of Ty as 



377(l + w B ) + 2 
277(1 + w B ) +4' 



(25) 



Therefore, for a nearly constant Ty, 77 and w<f,, the tracker 
condition, i.e., equation l|15|) gives the following possibil- 
ities: 

a. If < w Bl then increases with time. Then we 
conclude from equation (|15|l that \V ( f,/V^~ | decreases for 
a tracker solution. However, taking the time derivative 



of V+/V- 



d 
~dl 



we obtain 



V 2 

V+l 



VI 



V L 



Hence, {V^/V '^'l decreases if Ty 
w B is observed for 



> 



l+V 

2 ■ 



(26) 



Thus, w<f, < 



1 



V 



(27) 



Combining this with the condition for stable tracking be- 
havior fcauation 125(1 . we obtain 



Ty > max 



3t7(1 + w b ) + 2 I + 77 
277(l + w B ) + 4' 2 



(28) 



as it gives viable 
These conditions 



This is the most interesting case, 
models for an accelerating universe, 
encompass more solutions than the ones derived in Rcfs. 
[23L l24l |25| . For example, for the exponential potential, 
we have Ty — 1, and the above conditions are satisfied as 
long as 77 < 1 (including, for example, the Gauss-Bonnet 
expansion law). 

b. If w^ > w B , then tracking behavior is observed for 



3r](l + w B ) + 2 
277(1 + w B ) + 4 



< Ty < 



n- 



1 



(29) 



c. If Ty = (1 + t;)/2, then w$ = w B . This is one of the 
main results (using somewhat different notation) derived 
in Ref. [H. 



III. fc-ESSENCE 

A. Tracking solutions 

In general, fc-essence can be defined as any scalar field 
with non-canonical kinetic terms, but in practice such 
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models are usually taken to have a Lagrangian of the 
form: 

C = V{<t>)F{X), (30) 
where (f> is the scalar field, and X is defined by 

1. 



X = -V^V^ 



(31) 



The pressure in these models is given by 

P4> = C, (32) 

where C is given by equation l|3U|) , while the energy den- 
sity is 



V{<f>)[2XF x -F], 



(33) 



where Fx = dF/dX. Therefore, the equation of state 
parameter, = Pc/>/p<p, is just 



F 



W<f, 



2XF X - F 



(34) 



In defining the sound speed, we follow the convention of 
Garriga and Mukhanov [8j , who argued that the relevant 
quantity for the growth of density perturbations is 



Fx 



(dp/dX) 

(dp/dX) F x + 2XF XX ' 



(35) 



with F X x = d 2 F/dX 2 . 

In a flat Robertson- Walker metric, the equation of mo- 
tion for the fc-essence field takes the form: 

(F x + 2XFxx)4> + 3HF x <i>+(2XFx-F)¥t. = . (36) 

We can express the equation of motion for <fr in an alter- 
native form which will be useful for subsequent analysis: 



±^V2X = n( 1 -±p-) { u-, w, 



(37) 



where 



.4 : 



(XFx-F)(2XFxx+F x ) 
XFl - FF X - XFFxx 



l-w</, 



(38) 



y = (1 + w$)/(l — w^) and y' = d\i\y / d\na, and plus 
(minus) sign corresponds to </> < (<f) > 0), respectively. 
The tracker condition (w^ ~ constant) becomes 



± 



V A 



y(n+2)/2 



F 



i/2 



'2X 



(39) 



It is not surprising to see that the tracker condition for fc- 
essence has an extra "degree of freedom" in F(X). The 



functional form of F(X) plays a crucial role in deter- 
mining the tracking conditions for fc-essence and we shall 
consider it in the next section. 

After taking the time derivative of equation l|37|l and 
using equation (|12[) . we obtain 



2(1 + to*) 
[v(wb - w$)Q.b + V+(V - 2) w 4>]Ay' 



2(1 
2(1- 



- w<j,)(6 

w<p)y" 



Ay') 



{l + Wlj> )(6 + Ay>) 2 (c 2 -w^ 
2(w 4> (l-c 2 s )-{dc 2 s /dt)(l-w 4> ))y'/H 
{l+w^ + Ay') 2 (c 2 s -w^) 2 



(40) 



where y" = d 2 In y/d In a 2 . We note that for 77 = 1, equa- 
tion (|4*U|) reduces to the one derived in Ref. [l2j . 

For a background-dominated universe with a constant 
and almost constant Ty, the tracker equation l|40() 
reduces to 



n - 



2(1 + ^) : 

l + WB 



1 



Wtj, 



(41) 



Note that equation l|41|) for fc-essence closely resembles 
equation (|18|l for quintessence; the only difference is the 
constant appearing in the first term. For the standard 
Hubble expansion law (77 — 1), we obtain 



1 



1 + wb 



(42) 



l + w<p 

in agreement with the results of Ref. [l^ ■ 

B. Stability of the tracking solutions 

To determine the stability of the tracking solution, we 
repeat the calculation of Sec. II. B. for the case of fc- 
essence. We assume a fc-essence field with equation of 
state parameter wq and perturb wq by an amount d. 
Then we expand equation (|40|l to lowest order in 5 and 
its derivatives to obtain 



25 + 3(77(1 + w B ) - 2w }8 
+9r){l + w B )(c 2 -w Q )6 = 0, 



(43) 



where the prime means d/d\na. The solution of this 
equation is 



(44) 



where 



7 = -- [r/(l + Wb) - 2wq] 

± + ™b){c 2 - w ) - [77(1 + w B ) - 2w } 2 . 



(45) 
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Again, in order to have d decay, the real part of 7 has to 
be negative. Hence, it follows that 



w < 



77(1 + w B ) 



and 



w < c s 



(46) 



(47) 



At this point, the above conditions cannot be trans- 
lated into relations in terms of IV without considering 
the functional form of F(X), since w<f, and c 2 s both de- 
pend on F(X). Now we discuss the restrictions on the 
form of F(X) for constant w^. 

A variety of functional forms for F(X) and V{<f>) have 
been considered in fc-essence models (see, e.g., Refs. 
Ell ED) However, we will focus on the form of F{X) 
responsible for stable tracking solutions for a constant 
equation of state. In order to find the functional form of 
F{X) for stable tracking solutions with constant w^, we 
note that equation 1|34[1 can be written as 



dlnF(X) _ l + 
d\n X 2w$ 



(48) 



Case 1. The first possibility emerges if we treat equa- 
tion l|48|) as a differential equation and derive the general 
solution, which is 



F(X) = X?, 
where /3 is a constant, and is then 

1 



2/3-1 



(49) 



(50) 



By inserting equation (|49() into equation (|35|l . we obtain 



2/3 



so that 



(51) 



(52) 



These solutions were previously derived in Ref. |2jj ; we 
note here that they are independent of 77, and therefore 
of the expansion law. These solutions also do not depend 
on the form of V (</>). 

It is obvious from equation (|52|l that c 2 s < for any of 
these models with negative pressure (wj, < 0). If <? s < 0, 
then the fc-essence fluid is unstable against perturbation. 
Moreover, equation H49|) describes a phantom field for 
< /3 < 1/2. 

Case 2. A second class of solutions arises if the field 
evolves to a state for which X — Xq, where Xq is a 
constant 0,0. In this case, we have [!j 



d\nF(X) 



d In X 



1 



W0 



X=X 



2iU0 



(53) 



Again, we see that equation l|53|l is independent of 77 and 
hence, independent of the expansion law. However, the 
condition for a stable solution of the form X — Xq does 
depend on 77, as we now show. 

From equation l|41() . the tracking conditions, equations 
()46J) - H47JI . take the following form in terms of IV: 



and 



IV > 



IV > 1 



2n(l + w B ) + 2 
n(l + w B ) + 2 ' 

77(1 + 1115) 
2(1 + c*) ■ 



(54) 



(55) 



Therefore, for a nearly constant IV, 77, and w^, equation 
(I39|l gives the following possibilities: 

a. If < wb, then increases with time. Then 
we conclude from equation i|3l3|) that \\j2XV ( j > / F^V 2 ^ | 
decreases for a tracker solution. However, taking the time 

derivative of (y^xv^/Fiv^ , we obtain 



'2X 



dt \ v 1 



Fi 



2XV[ 

Fi v v -i 



(56) 



In the derivation of this equation, we have used the 
condition that X = X . Hence, | ^IXV^j F^V^ \ de- 
creases if T v > (rj + 2)/2. Thus, w,p < wb for 



V 



(57) 



Combining this with the conditions for stable tracking 
behavior ( equations I54I55|I . we obtain 

r v >ma X [— T){1+WB) + 2 ^ + W ^y}- (58) 
b. If ?«0 > wb, then tracking behavior is observed for 



[ 2r ] (l + w B ) + 2 g(l + wg) 

max\ ; r , 1 H ; -— 

L r){l + w B ) + 2' 2(1 + C 2) J 



< Ty < 



V 



(59) 



c. If Ty = (r/ + 2)/ 2, then = wb- This case en- 
compasses the solutions presented in Ref. pB| . 



IV. DISCUSSION 

We have extended the formalism in Refs. @ and 0] 
to derive the tracker conditions for quintessence and fc- 
essence, respectively, for an arbitrary cosmological ex- 
pansion law, H 2 = f(p), when the universe is dominated 
by a background fluid. Our main new result is that, with 
the exception of the special cases discussed above, track- 
ing solutions for either quintessence or fc-essence are pos- 
sible only for r] = d\nf /d\np w constant, which is the 
case only when /(p) is well-approximated as a power-law. 
In fact, such power-law behavior corresponds to most of 
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the models previously considered for non-standard ex- 
pansion laws. 

We note further that the expressions for w$ for both 
quintessence and fc-essence, and the conditions for stable 
tracking behavior, can be derived by replacing 1 + wb 
with 77(1 + wb) in all of the corresponding equations for 
the standard expansion law. This is not surprising, since 
a given value of wb corresponds to a background density 
scaling as ps oc a~ 3l - 1+WB \ Taking a constant value of 77 
in equation 10 then gives H 2 oc a~ 3v ( 1+WB \ so 1 + wb is 
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